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ABSTRACT 

Most of the physical problems including sound waves in a viscous medium, waves in fluid filled viscous elastic 
tubes and magneto hydrodynamic waves in a medium with finite electrical conductivity are modeled by nonlinear partial 
differential equations. Many numerical and analytical methods are used to solve non-linear partial differential equations. 
Wavelets have generated a huge interest in different areas of applied mathematics, physics and engineering. 
Wavelets have been applied to the numerical solution of partial differential equations. In this paper we have developed a 
wavelet Galerkin scheme for non-linear partial differential equations. The present scheme is efficient, accurate and has got 
several advantages over other numerical methods. 
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INTRODUCTION 

Wavelet theory [4,5] was originally applied as a powerful tool for signal and image processing. 
Wavelet analysis allows to obtain an efficient sparse representation of some function particularly good localization both in 
space and in scale of each element of a given wavelet basis. In fact, in the wavelet expansion of a function many 
coefficients are negligible and by discarding these coefficients we can obtain a sparse but accurate approximate 
representation. Moreover, using wavelet decomposition it is possible to detect singularities, irregular structure and transient 
phenomena exhibit by the analyzed function. In the last years, the good features of wavelets have generated a huge interest 
in different areas of applied mathematics, physics and engineering. Wavelets have been applied to the numerical solution 
of partial differential equations (PDEs). These kinds of methods have been studied both from the theoretical and the 
computational point of view. One of the main topics arising in the application of wavelets to the numerical algorithms for 
PDEs is the study of adaptivity solutions to PDEs often behave differently in different areas of their domain [4,5]. 

WAVELETS 

A wavelet [16] is a wave-like oscillation with amplitude that starts out at zero, increases, and then decreases back 
to zero. It can typically be visualized as a "brief oscillation" like one might see recorded by a seismograph or heart monitor. 
Generally, wavelets are purposefully crafted to have specific properties that make them useful for signal processing. 
Wavelets can be combined, using a "shift, multiply and sum" technique called convolution, with portions of an unknown 
signal to extract information from the unknown signal. 
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TYPES OF DISCRETE WAVELET TRANSFORM 

The following are the different types of discrete wavelet transform [9] 

• Haar Wavelet 

• Beykin 

• BNC wavelets 

• Coiflet 

• Cohen-Daubechies-Feauvean 

• Daubechies Wavelet 

• Binomial-QMF 

• Legendre Wavelet 

• Villasenor Wavelet 

• Symlet 

DAUBECHIES WAVELET 

Daubechies wavelets [7,9], based on the work of Ingrid Daubechies, are a family of orthogonal wavelets defining 
a discrete wavelet transform and characterized by a maximal number of vanishing moments for some given support. 
With each wavelet type of this class, there is a scaling function which generates an orthogonal multiresoution analysis. 

Daubechies wavelets are chosen to have the highest number of A of vanishing moments for given support width 

N=2A and among the 2 A 1 possible solutions one is chosen whose scaling filter has extremal phase. Daubechies wavelets 
are widely used in solving a broad range of a problems, e.g. self- similarity properties of a signal or fractal problems, signal 
discontinuities, etc. Daubechies orthogonal wavelets D1 — D20 (even index number only) are commonly used. 
The index number refers to the number ./V of coefficients. Each wavelet has a number of zero moments or vanishing 
moments equal to half the number of coefficients. A vanishing moment limits the wavelet's ability to represent polynomial 
behavior or information in a signal. 

HAAR WAVELET 

In mathematics, the Haar wavelet [9] is a sequence of rescaled "square-shaped" functions which together form a 
wavelet family or basis. Wavelet analysis is similar to Fourier Analysis in that it allows a target function over an interval to 
be represented in terms of an orthonormal function basis. The Haar sequence is now recognised as the first known wavelet 
basis and extensively used as a teaching example. 

The Haar sequence was proposed in 1909 by Alfred Haar . Haar used these functions to give an example of an 
orthonormal system for the space of square-integrable functions on the unit interval [0, 1]. The study of wavelets, and even 
the term "wavelet", did not come until much later. As a special case of the Daubechies wavelet , the Haar wavelet is also 
known as D2. The Haar wavelet is also the simplest possible wavelet. The technical disadvantage of the Haar wavelet is 
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that it is not continuous , and therefore not differentiable . This property can, however, be an advantage for the analysis of 
signals with sudden transitions, such as monitoring of tool failure in machines. 

The haar wavelet's mother wavelet function l/f(t) can be described as 



1 0<f<l/2, 
-1 \/2<t<\, 
0 otherwise 



Its scaling function (j) ( t ) can be described as 



|1 0< t <l, 
1 0 otherwise. 



WAVELET GALERKIN METHOD 

The Galerkin method belongs to the family of weighted residuals methods, where the solution to a partial 
differential equation, u [x, t ) is approximated by a finite series of functions (J) k ( x) as follows: 

p 

« O'0=X a k (Oh ( x ) 

(1) 

where <f> k (x)are the basis or trial function, a k (?)are coefficients to be determined (possibly time-dependent) 
that satisfy the partial differential equation, and p are the number of functions. In general the approximate solution does 

not satisfy the partial differential equation exactly, and substituting its value results in a residual, which in turn is 
minimized in some sense. 

Consider the following boundary value problem (BVP)[12] 
Lu(x,t)- f (x,t) = 0 ,xeQ.,t>0 

(/■) 

Substituting the trial solution (1), we get 

p 

L E fl *(0&(*)-/U0=*(*.0 

*=i (3) 
The method of weighted residuals minimizes R[x, ?)by forcing it to be zero in the domain £1 using weight 
functions w ■ ( x) such that, for every weight function, 

^ R(x,t)wj(x)d x=0 for 7=1,2, ,q 

n 

where q is the number of weight functions to be determined, depending on the boundary conditions and number 
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of scaling functions. This discretization leads to a system of linear or non-linear ordinary differential equations where the 
values for a k (?) can be determined. In the Galerkin method, the weight functions are chosen to be the basis functions of 
Daubecheis wavelet. 

In wavelet Galerkin method we choose wavelet basis as weight functions. Having multiresoltuion analysis, 
Vj , j €E Z with scaling function (f) (jc) , one can use <f> . k ( x) as the basis functions for the Galerkin method. We know that 



one can 



the set forms an orthonormal basis of Vj . Thus at approximation level j 

take the orthogonal projection of u onto in the following manner 

P . . 

u(x,t)**r.u(x,t) =2X* (0^,* (4 a j,k (0 =\ u M^m M) 

*=i 1 (4 ) 

and this will provide an accurate approximation to u . Furthermore, for some j , V, will capture all details of the 
original function. 

WAVELET GALERKIN SCHEME 

PROBLEM 1: Consider the modified Burgers 'equation (MBE)[13] is based upon Burgers' equation of the form 

U t + U 2 U X - VU XX - 0 (5) 

with initial condition 

u(x,0) = f{x) (6) 
and boundary conditions 

u x (0,t)-u x (l,t)=0 

««(o,0= M «(i.0=° ( ? ) 

The equation has the strong nonlinear aspects and has been used in many practical transport problems, for 
instance, nonlinear waves in a medium with low-frequency pumping or absorption, turbulence transport, wave processes in 
thermo-elastic medium, transport and dispersion of pollutants in rivers and sediment transport, and ion reflection at 
quasi-perpendicular shocks. 

First we develop a weak formulation, from which we will derive the discretization. Multiplying both sides of 
Burger's equation by a test function, (j). e V and integrating, we get 

1 1 1 

\ u , 0j, P W d x + J u 2 u x 0 j p (x)dx-v^u xx 0 jtp (x)dx = 0 (8) 
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Let the trial solution of partial differential equation (8) be 

L 

u(x,t) = Y j a k {t)(j) jk (x) (9 ) 

where <j> . t (x) = 2^ </>(2 j X-k) and 0(x) is a scaling function of Daubechies compactly supported wavelet 

with vanishing moment N j 2 . Now taking test function as ^ (x) = 2^ 2 fy[2J X— /?) , we get a system of simultaneous 
differential equation. 

1 L 1 L 



{Z«* ( x M> (*)^+4Z a * ( x )IX (*M> (*)Z a « (OC ( x H> 

o >t=o o *=o m=a «=o 

1 L 

+ w J Z a * ( r ) ( x ) ( x ) = 0 



k=Q 



L L L ' 

^ (0 =-Z fl * (0Z a <» (0Z a « (OjA* ( x M» ( x )C ( x M> ( x ) dx 

fc=0 m=0 re=Q o 



*=0 o 



% (t)=-t * (0 Z (OZ - *>Z «* ('K; 1 , do) 



fc=0 m=0 n=0 k=0 

where, 

l 



"KS, = JV„* ( X M> ( X )C ( X M„ ( X )^ X C») 

0 

1 

(12) 

0 

PROBLEM 2: Consider Kortewege-de Vries-Burgers' equation [14] is a nonlinear partial differential equation, 
which is given by 

U t + SUU x - V« xv + jUU xxx =0 (13) 
with the initial condition 

u(x,0) = u Q (x) (14) 
and boundary conditions 
U x (0,t) = U x (l,t) = 0 
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"^(o>0 =M ^( 1 '0 =0 ( 15 > 

where £, v and ju are positive parameters. This equation was initially formulated by Gardner [2]. 
This model arises in many physical applications such as propagation of undular bores in shallow water waves, propagation 
of waves in elastic tube filled with a viscous fluid and weakly nonlinear plasma waves with certain dissipative effects. 

It represents long wavelength approximations where effects of the nonlinear advection term UU X are counterbalanced by 

the dispersion U . A number of theoretical issues related the KdVB equation have received considerable attention. 

The traveling wave solution in particular has been studied extensively. Demiray [6], Antar and Derniray [11] derived 
KdVB equation as the governing evolution equation for wave propagation in fluid-filled elastic or viscoelastic tubes in 
which the effects of dispersion, dissipation and nonlinearity were present. Eq.(13) is combination of the Burgers' equation 
{/I = 0 ) and the KdV equation (V = 0). Burgers' equation [2] was first used by Burger for the study of turbulence in 

1939, whereas KdV equation [2] was first suggested by Kortewege and de Vries who used this as a nonlinear model to 
study the change in shape of long waves moving in a rectangular channel. KdVB equation has been solved by many 
authors exactly and numerically. 

We develop a weak formulation, from which we will derive the discretization, multiplying both side of KdV 
Burgers equation by a test function, 0. ( x} £ V and integrating we 

11 11 

J u,<f>j, P {x)dx + e\ uuj> jp (x)dx- vj uj> jp (x) dx+/i\ u m <p jp (x)dx = 0 (16) 

0 0 0 0 

Let the trial solution of partial differential equation (13) be 

L 

u(x,t) = Y j a k {t)(j) jk {x) (17) 

k=Q 

where 0 . k {x) = 2/ 2 X~k) and 0(x) is a scaling function of Daubechies compactly supported wavelet 

with vanishing moment N/2. Now taking test function as^ . [x) = 2 <I)(7J X — we get a system of simultaneous 
differential equation. 

1 L 1 L 



1^ 1^ 

{t)0j, k {x)0j, P {x)dx + e\Y J a k MIX (0#.m i x )hp ( x ) dx 

0 i=0 o A '=° m=Q 

1 L 1 L 

+v\^a k (f)^ {x)4> hp {x)dx+{i]Y,a k (t)</)' jk {x)^ p {x)dx = 0 

0 *=0 o 
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L L 1 

k=0 m=0 o 

L 1 L 1 

(OJVm ( ? )JX* ( x ) dx 



k=0 o o 



s (0 - -*I>* (05>. ('Ki (0 a 1'!p -/<I>* as) 



k=0 m=0 lfc=0 A=0 

where, 



^Z=^jA x )^ m ( x )h P ( x ) dx < 19 > 

0 

1 

< 20 > 

0 
0 

PROBLEM 3: Consider the Newell- Whitehead- Segel equation [13] models the interaction of the effect of the 
diffusion term with the nonlinear effect of the reaction of the reaction term. The Newell-Whitehead- Segel equation [1] is 
written as: 



du d Z 



U +u[\-u 2 ) (22) 



dt dx 2 
with the initial condition 

u(x,0) = u 0 (x) (23) 
and boundary conditions 

U x (0,t) = U x (l,t) = 0 (24) 
Let the trial solution of the partial differential equation (22) be 

L 

K(x,0=Z fl *('M,*(*) (25) 

k=0 

Now integrating both sides of equation (22) with respect to weight function (J), (jc)we get a system of 
simultaneous differential equations of first order 

ju t vdx- ju xx vdx + juvdx- juu 2 vdx 
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Jw,^ p (x)dx = Jw vv ^. p (x)dx + §u0j ■ p (x)dx -^uu 2 (f) ] p (x)dx 

L L 

\Y. a 'k ( f )hk ( x )0j,p ( x ) dx = \H a k (0^1.* ( x )0j.p ( x ) dx 

k=0 k=Q 

+j"Z fl * (*M-,„ {Z a * (0 fl „, {*)<*, (0^/,* (*M> ( x ) dx 

s (0 = I>* (OA, 2 : 0 , - 1 (*K (0«. (Onffii +«, (0 < 26 > 



where 



Aj[J = J (27) 



and 

"SSi = R* (*M>, (*)*,., (28) 

-/ 

CONCLUSIONS 

In the present paper a new wavelet Galerkin scheme is developed for solving non-linear PDEs. 
Although a Wavelet Galerkin scheme is shown using Daubecheis wavelet, the same can be extended by using Haar 
Wavelet and coiflets Wavelets. The scheme is shown for three well-known non-linear PDEs which are arising in 
engineering and science. However it can be easily extended to other types of nonlinear PDEs. In the traditional finite 
element method the approximating functions are chosen in terms of polynomial over a finite element. With the large 
number of elements, the computational cost will be higher. To overcome this difficulty, the approximating functions are 
replaced with wavelet functions which reduce the computational cost. The present scheme is effective, efficient and can 
handle large number of PDEs. 
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